Evidence of Charge Density Wave Ordering in Half Filled High Landau Levels 
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We report on numerical studies of two-dimensional electron systems in the presence of a perpen- 
dicular magnetic field, with a high Landau level (index N > 2) half filled by electrons. Strong 
and sharp peaks are found in the wave vector dependence of both static density susceptibility and 
equal-time density- density correlation function, in finite-size systems with up to twelve electrons. 
Qualitatively different from partially filled lowest (JV = 0) Landau level, these results are suggestive 
of a tendency toward charge density wave ordering in these systems. The ordering wave vector is 
found to decrease with increasing N. 
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Two-dimensional (2D) electron gas systems subject to 
a perpendicular magnetic field display remarkable phe- 
nomena, reflecting the importance of electronic correla- 
tions. The most important among them is the fractional 
quantum Hall effect (FQHE), which was found in the 
strong field limit, where the electrons are confined to the 
lowest (N = 0) or the second (N — 1) Landau levels. The 
physics of FQHE is reasonably well understood p[ : the 
kinetic energy of the electrons are quenched by the strong 
perpendicular magnetic field and the Coulomb interac- 
tion dominates the physics of the partially filled Landau 
level; at certain Landau level filling factors [y, defined to 
be the ratio of the number of electrons to the number of 
Landau orbitals in each Landau level) the electrons con- 
dense into a highly-correlated, incompressible quantum 
fluid, giving rise to quantized Hall resistivity (p X y) and 
thermally activated longitudinal resistivity (p xx ). 

Experimentally, the FQHE has never been found at 
filling factors v > 4, when the partially filled Landau 
level has Landau level index N > 2 (taking into account 
the two spin species of the electrons). Nevertheless, re- 
cent experiments |^,^| on high quality samples have re- 
vealed remarkable transport anomalies for v > 4, espe- 
cially when v is near a half integer, which means the 
partially filled Landau level is nearly half filled. Such 
anomalies include a strong anisotropy and nonlinearity in 
p xx . They reflect intriguing correlation physics at work 
in these systems that is qualitatively different from the 
FQHE and yet to be completely understood. 

It was argued before the discovery of the FQHE, 
that the ground state of a 2D electron gas in a strong 
magnetic field may possess charge density wave (CDW) 
order. Recent Hartree-Fock (HF) calculations ||-|7j find 
that single-Slater determinant states with CDW order 
have energies lower than the Laughlin-type liquid states 
for N > 2. The CDW state with ID stripe order, or 
stripe phase J5|-fj]], which is predicted to be stable near 
half filling for the partially filled Landau level, can in 
principle give rise to transport anisotropy as the orienta- 
tion of the stripe picks out a special direction in space. 



Questions remain, however, with regard to the stability 
of the HF states against quantum fluctuations as well as 
disorder, especially when N is not too large. 

In this paper we report on results of numerical studies 
of partially filled Landau levels with N > 2, in finite size 
systems with up to N e = 12 electrons and torus geometry. 
The torus is ideal for this study, because, while it main- 
tains translational and rotational invariance in the plane 
for the infinite system, it does break them for finite sizes 
and therefore produces a preferred direction that allows 
for the CDW state to be aligned in one direction. The 
spherical geometry (which is another popular geometry 
for finite size studies), on the other hand, will rotationally 
average the state and necessarily introduce defects. We 
assume the magnetic field is sufficiently strong so that 
the filled Landau levels are completely inert and mixing 
between different Landau levels can be neglected. We 
calculate numerically the energy spectra, the wave vec- 
tor dependence of the static density susceptibility (x(l)) 
and the density-density correlation in the ground state 
(So (<?)). We find strong and sharp peaks in both x(<?) 
and 5*0 (g). These results are strongly suggestive of a ten- 
dency toward charge density wave ordering in the ground 
state. In accordance with this, we find nearly degener- 
ate low-energy states that are separated by the ordering 
wave vector. 

The methods used in the present study are identi- 
cal to those used in numerical studies of the partially 
filled lowest Landau level J8|. Since the kinetic energy is 
quenched by the magnetic field, the Hamiltonian contains 
the Coulomb interaction alone, which, after projecting 
onto the Nth Landau level |^], takes the form 

ff = ^^ e -« 3 /2 [jLjv( ^ /2)]2y(g)eiq .(R i -R 3 .) ; (1) 
i<j q 

where is the guiding center § coordinate of the 
ith electron, L^{x) are the Laguerre polynomials, and 
V(q) = 2ire 2 /q is the Fourier transform of the Coulomb 
interaction. The magnetic length I is set to be 1 for 
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convenience. In this work we impose periodic boundary 
conditions in both x and y directions in the finite-size 
system under study, and the q's are wave vectors that 
are compatible with the size and geometry of the sys- 
tem. We diagonalize the above Hamiltonian exactly and 
calculate various correlation functions. 
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FIG. 1. Energy (in unit of e 2 /el) versus the x-component 
of the momentum for the half-filled N = 2 Landau level with 
ten electrons, rectangular geometry, and aspect ratio 0.75. 
The momenta of the five nearly degenerate low-energy states 
are (±0.485, 0), (±1.456, 0), and (2.427, 0). 

In Fig. 1 we present the energy spectrum of half-filled 
N = 2 Landau level with ten electrons, rectangular geom- 
etry, and aspect ratio 0.75. The spectrum is qualitatively 
different from the known incompressible FQHE states in 
the following aspects: i) the ground state is not separated 
from the excited states by a gap; instead, there are several 
nearly degenerate low energy states, separated by a char- 
acteristic wave vector q* (the physical importance of q* 
will be discussed later). This almost exact degeneracy is 
not specific to this particular geometry; in Fig. 2 we have 
plotted the energy levels (for all momenta) versus aspect 
ratio a, in which it is clear that the near degeneracy is 
present for 0.7 < a < 1.0. A gap in the spectrum, on the 
other hand, is the most important property of a FQHE 
state that gives rise to incompressibility. ii) The momen- 
tum of the ground state is, in general, not related to any 
reciprocal lattice vectors and is sensitive to the geometry 
of the system; such sensitivity to boundary conditions is 
characteristic of compressible states. The known FQHE 
states, however, all have ground state momentum equal 
to one half of a reciprocal lattice vector and are indepen- 
dent of system geometry, reflecting intrinsic topological 
properties of the state. The spectra and quantum num- 
bers of these ground states are also very different from the 
Fermi-liquid like compressible state at half filling in the 
lowest Landau level 101. In that case cluster- like ansatz 



variational wave functions pfXflll have remarkably large 
overlaps with the exact ground state and, from the shape 
of the optimal cluster for a given geometry, predict the 
momentum of the ground state with remarkable accuracy 
[fyj. The momenta of the ground states here, however, 
do not match those of the optimal cluster states. Impor- 
tant differences are also seen in the response functions 
studied below. We thus conclude that these states are 
compressible and have properties that are qualitatively 
different from the Fermi-liquid like compressible state at 
half filling in the lowest Landau level. 
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FIG. 2. Energy levels versus aspect ratio for half-filled 
N = 2 Landau level with ten electrons and rectangular ge- 
ometry. The inset is a blow-up of the low-energy spectra for 
aspect ratio between 0.7 and 1.0; the energies have been mul- 
tiplied by 10 4 . 

We now turn to the discussion of response functions. 
Here the fundamental quantity of interest is the dynam- 
ical structure factor, defined to be 

S (q,w) = ^ £ |(0| ^e^|n)| 2 ^„ -E -co). 



(2) 



The summation is over states in the given Landau level 
that is being studied. Various physically important quan- 
tities may be calculated from So(q, u>) Q; in particular, 
the projected static density response function is the in- 
verse moment of 5rj(q, UJ ) : 



x(q) 



dujSo(q, uf)/oj, 



(3) 



and the projected equal time density-density correlation 
function is the 0th moment of So(q, w): 



So(q) 



dw5 (q,w). 



(4) 
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FIG. 3. So(q) and X (q) (in unit of el/e 2 ) for half-filled 
N = 2 Landau level with ten electrons. The data were taken 
from systems with rectangular geometry, and a number of 
different aspect ratios. The highest peak value for 5o(q) and 
x(q) are 3.7 and 39726.3 respectively, which occur at aspect 
ratio 0.75. 

In Figs. 3 and 4 we present the q dependence of x(q) 
and So(q) for half filled N = 2 and N = 3 Landau levels 
respectively. The data were taken from systems with 10 
electrons and rectangular geometry, for several different 
aspect ratios. We see sharp and strong peaks in x(q) at 
q* = (0.97, 0) and q* = (0.84, 0) for N = 2 and N = 3 
Landau levels respectively. At the scale of the peak value 
of Xi t ne response at different q's are indistinguishable 
from zero, except for a secondary peak positioned at ex- 
actly 3q*, the height of which is much smaller but quite 
noticeable. 

This is strongly suggestive of a tendency toward CDW 
ordering, as an (unpinned) CDW responds strongly to an 
external potential modulation with a wave vector that 
matches one of its reciprocal lattice vectors. The fact 
that secondary peaks appear only at integer multiples 
of the primary wave vector (q* ) suggests that the CDW 
has ID (stripe) structure. The absence of response at 
2q* (and any other even multiples) is consistent with the 
presence of particle-hole (PH) symmetry in the underly- 
ing Hamiltonian at half filling: the PH transformation 
of the CDW state is equivalent to translation by half a 
period. Consistent with the stripe CDW picture, sharp 
peaks are also observed in So(q) at q = q*, indicating 
strong density-density correlation in the ground state at 
the ordering wave vector. Similar behavior is also seen 
at higher Landau levels and other filling factors. Again, 
these results are in sharp contrast with FQHE states (like 
v = 1/3 in the lowest Landau level), or the Fermi liquid- 
like state of the half filled lowest Landau level. In the 
former case, x is small (below 20) for all g's due to the 



large gap in the spectrum || ; while in the latter case al- 
though there are peaks in both x(o) and So(q), the peak 
position is tied to q — 2kp — 2 and the height of the 
peak in x( a ) is i n general much smaller than those seen 
here and the peaks are broader pT| . 




FIG. 4. Same as in Fig. 3, for N = 3 Landau level. The 
highest peak values are 4.17 and 2953432.4 respectively, which 
occur at aspect ratio 0.56. 

The origin of the strong peaks in x at q* may be traced 
to the almost exact degeneracy of the low energy states in 
the spectra that are separated by q*. Since these states 
are connected by a potential modulation with wave vector 
q* , an extremely small energy denominator ensures huge 
response from the system. We emphasize that such near 
degeneracy is a generic property of the system and not 
specific to a particular system size or geometry Jl2| . This 
can be seen clearly in Fig. lb, where we plot energy levels 
from different geometries together. If the system devel- 
ops long-range CDW order in the thermodynamic limit, 
we would expect the number of nearly degenerate states 
to increase with system size; in the thermodynamic limit, 
there will be infinitely many of them, spaced in momen- 
tum by q* , that become exactly degenerate. In this case a 
ground state that spontaneously breaks the translational 
symmetry may be constructed by taking linear combina- 
tions of these degenerate states, even though this is not 
possible (unless degeneracy is exact) in any finite system 
where all eigenstates must have a definite momentum. 

In Fig. 5 we plot the Landau level dependence of q* for 
half filled Landau levels. For systems with a given num- 
ber of electrons in a given Landau level, we define q* to 
be the modulus of the wave vector that gives the largest 
X(q) for all geometries (aspect ratios). Despite the small 
and non-systematic dependence on the number of elec- 
trons, it is clear that q* decreases with increasing N. We 
are unable to accurately determine q* beyond N = 4, as 
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in this case 1/q* becomes comparable to (and eventually 
exceeds) the linear size of the largest size system that we 
are able to study. This result is qualitatively consistent 
with the prediction of Hartree-Fock theory |5[j^], which 
predicts that the period of the CDW is set by the scale 
of the cyclotron radius in a given Landau level. In units 
of inverse magnetic length, this implies q* oc 1/ y/2N + 1. 
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FIG. 5. Estimated ordering wave vector q* versus Lan- 
dau level index N. The solid curve is the prediction of 
Hartree-Fock theory. 

We have also performed numerical studies away from 
half-filling, at filling factors v = 1/4,1/3,2/5, etc, for 
Landau levels 2 < N < 6. At all these filling factors, 
we have seen qualitatively similar behavior in \ an d So, 
which are suggestive of a tendency toward CDW order- 
ing. No evidence of incompressible FQHE states is found. 
This is consistent with Hartree-Fock theory, which pre- 
dicts H| no FQHE for N > 2. Hartree-Fock theory 
also predicts that when sufficiently far away from half- 
filling, the stripe phase becomes unstable against a differ- 
ent type of CDW ordering, the "bubble phase" fl , which 
has a two-dimensional lattice structure. In our calcula- 
tions we have seen some indication that this may be the 
case; however more work is needed to make a clear dis- 
tinction between these two types of structures. We leave 
this to future investigation. 

As stated above, the physical properties of the sys- 
tems that we have studied here are qualitatively consis- 
tent with predictions of Hartree-Fock theory. We also 
find that in the single Slater determinant basis for the 
wave functions, the Hartree-Fock single Slater determi- 
nant (with simple stripe structure) has the highest weight 
in the exact ground state. For example, in N = 2 and 
N = 3 Landau levels with ten electrons and rectangular 
aspect ratios of 0.75 and 0.56 respectively, the highest 
weight single Slater determinant has the occupation pat- 



tern (in Landau gauge) 1 1 1 1 1 00000 1 1 1 1 100000 (and ones 
that may be obtained by translating this pattern) , where 
1 stands for an occupied orbital and stands for an empty 
orbital. The maximum weights are 0.1463 and 0.1986 for 
N = 2 and N = 3 Landau levels; the next highest weight 
single Slater determinant has a weight that is approxi- 
mately fifteen times smaller for N = 2 and 1400 times 
smaller for A = 3. However, by making a linear combina- 
tion of the five nearly degenerate states we can single out 
the above occupation pattern and construct an approxi- 
mate eigenstate that breaks translational symmetry. The 
overlaps (squared) with HF wave function then become 
0.7338 for N = 2 and 0.9930 for A = 3. While there are 
still some fluctuations on top of HF states for N = 2, 
these fluctuations are completely gone for A = 3. Frad- 
kin and Kivelson IS] recently considered the effects of 
thermal and quantum fluctuations on the Hartree-Fock 
stripe phase, and predicted a number of novel phases. 
Due to the limited system sizes in numerical studies, we 
are unable to distinguish among these phases (whose dis- 
tinctions show up at large distances only). 
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